Exact solutions in bouncing cosmology 
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We discuss the effects of a (possibly) negative (1 + z) 6 type contribution to the Friedmann equation 
in a spatially flat universe. No definite answer can be given as to the presence and magnitude of a 
particular mechanism, because any test using the general relation H(z) is able to estimate only the 
total of all sources of such a term. That is why we describe four possibilities: 1) geometric effects of 
loop quantum cosmology, 2) braneworld cosmology, 3) metric-affme gravity, and 4) cosmology with 
spinning fluid. We find the exact solutions for the models with p 2 correction in terms of elementary 
functions, and show all evolutional paths on their phase plane. Instead of the initial singularity, the 
^ ' generic feature is now a bounce. 

o ' 

_ 1 I. INTRODUCTION 

pLn \ Cosmology, mainly due to its modern astronomical observations, seems to have been limited to what is nowadays 
dubbed the Cosmological Concordance Model (CCM). In the observational cosmology this role is played by the cold 
dark matter model with the cosmological constant (ACDM model). It offers the simplest explanation of the current 
(f") . Universe filled with two components-dark matter and dark energy [l], Q ■ The former is in the form of non-relativistic 
dust matter contributing one third of the total energy density. The latter has negative pressure, violating the strong 
energy condition, and constitutes the remaining two thirds. 

Interpretations of observational data within this model lead to the conclusion that we live in an accelerating universe 
which is almost spatially flat and low density 0, 0,0. Although the introduction of A could be justified on purely 
geometrical ground, many authors seek a physical source of a suitably behaving matter, mimicking such a geometrical 
effect. In literature there are many cosmological models with dark energy which take part in the contest for the best 
. description of the accelerating Universe @, 0] • 

On the other hand, the very early universe, when classical physics is no longer sufficient, is not so well understood. 
^jH, It is believed that "initial conditions" for the classical evolution could be obtained from different, fundamental physics 
which dominates when the energy density is high enough, and which turn negligible as the Universe expands. In other 
j^Tjl words, new physics is welcome in the process of cosmological model building [H, [^j . In particular, one could hope for 



avoiding the initial singularity with a bounce of loop quantum cosmology [1C 

What we consider here are the effects of theories which predict a p 2 -type modification to the Friedmann equation, 
which fits the above description. Although such reduction is a big simplification, we feel it is important to stress that 
such problems are explicitly solvable, which is often not checked before applying numerical studies in current works; 
and that it is also possible to find and classify all types of evolutional scenarios which help to gain insight into more 
complicated models. 

Below, we outline some recently significant possibilities of non-classical physics, and proceed to solve the resulting 
Friedmann equation in the subsequent sections. 



A. Braneworlds 



The basic idea is that our observational universe is some four dimensional surface (called a brane) embedded in 
a more dimensional bulk spacetime, in which the gravitational field, but not the others, can freely propagate. The 
universe is self accelerating due to an additional term appearing in the Friedmann equation when constrained to the 
brane (for a pedagogical introduction into extra dimensions cosmology see e.g. [Til 

Let us consider a higher dimensional cosmology in the Randall-Sundrum framework 
equation assumes the following form [ll| [ljj 

where A4 is the four dimensional cosmological constant, e = ±1, and c is an integration constant whose magnitude as 
well as sign depend on the initial conditions. The fifth term (called dark radiation if c < 0) we put equal to because 



14, njl . Then the Friedmann 
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the third term is of the main interest in the present discussion. This p 2 contribution arises from the imposition of a 
junction condition for the scale factor of the brane. Note that the p 2 term would decay as rapidly as a -6 in a matter 
dominated universe, thus modifying the early evolution and not being significant later on. 

Both negative and positive e are possible because e corresponds to the metric signature of the extra dimension [l6| . 
However, the sign of e is crucial because models with timelike extra dimension can avoid initial singularity by the so 
called "bounce" [HI • Of course, the sign of e remains an open question - in [l9| , the authors discuss some consequences 
of the choice e = — 1. Then, the presence of the p 2 term on the right-hand side of ([T|) leads a contracting universe to 
a bounce instead of a big bang type of curvature singularity (p — > oo, R a bcdR abcd —* oo). If that happens, we have 
a constraint on the value of brane tension \a\ > (lMeV) 4 because the bounce takes place at densities greater than 
during nucleosynthesis. 

It should also be added that string theory does not provide any reasons for an additional timelike dimension. It is, 
nevertheless, a mathematically viable possibility, that some authors choose to investigate. 

Note that such modification was also recently investigated in the context of two-brane model [20( , and such equations 
arise in classically constrained gravity of Gabadadze and Shang [HI, [HI • 



B. Loop quantum universes 



The standard Big Bang cosmology presents us with the problem of initial conditions for the Universe. What we can 
observe from "the inside" cannot answer the question why particular conditions were chosen, as we are simply looking 
at the dynamics. Unfortunately we cannot look from "the outside" and risk any kind of statistical or other external 
explanation. Many authors shift the problem to the Planckian epoch, during which the quantum effects are important. 
A contraction phase preceding the present expansion, as predicted by some of these theories, seems to be one of the 
more attractive scenarios. It is very interesting that the geometric effects in Loop Quantum Cosmology predict a — p 2 
modification to the Friedmann equation [23l . |24| . This modification is relevant in the very early universe when its 
matter energy density becomes comparable with the Planck density. The effective Friedmann equation becomes 

, 8nGp ( p \ A 
H 2 = ^ 1-— +77. ( 2 ) 

3 V Pent/ 3 

where p cr it = v / 3/(167T 2 '7 2 )ppi HI, and one can distinguish classical and quantum bounces depending on the relative 
magnitude of p and p C rit- 

This is an application of the LQG methods, where the spacetime is discrete on the quantum level, considerably af- 
fecting the large scale [2(| [27|, Hll with nonperturbative corrections. The p 2 modification is shared with the braneworld 
model, and some comparison work can be found in see also 13211. 



C. Non-Riemannian cosmologies 

This misleadingly dubbed concept does not in fact involve the abandonment of Riemannian manifold structure, but 
the name has been widely used so far and we also adopt it here. 

The first non-Riemannian cosmological models were based on the Einstein-Cartan theory, which is a modification 
of general relativity by adding the torsion of the space-time 33]. In such models, the effects of spin and torsion 
are manifested by the presence of an additional — p 2 oc — a~ 6 contribution in the Friedmann equation. The main 
motivation for studying such models was the problem of the initial singularity which could be avoided due to the spin 
effects. 

Among various extensions of general relativity, the so called metric-affine gravity (MAG) is recently of interest. In 
contrast to the Riemann-Cartan theory, the connection is no longer metric which implies that the covariant derivative 
of the metric does not vanish. For a review see e.g. [1, |34j. 

This time, the Friedmann equation reads 

2 8nGp A ^ 2 
^=^+3+^' (3) 

where the new constant v can be of both signs, and ip is an integration constant. Therefore, the non-Riemannian 
quantities of this model (torsion and non-metricity) will modify very early stages of evolution and are negligible at 
later times. 

We find different possible interpretations of the presence of the p 2 contribution in the Friedmann equation. However, 
we must remember that cosmography which is based on the behaviour of the null geodesies maps the geometry and 
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kinematics of the universe in terms of H{z), without any reference to the source of each contribution. This happens 
because it measures only the average properties of the matter density, and eventually only the overall term (— )(l + z) 6 
is significant [35j ] - 

Interestingly, similar a~ 6 modifications have also been obtained in universes with varying constants (36j . It was 
shown that in the generic case, alpha varying models lead to a bouncing universe. 



II. DYNAMICS OF THE MODEL 



We take the Friedmann equation in the following form 

a'{tf 8ttG , n2 A 

^F = — M +V M + 3< (4) 

with a general p 2 term, and a new, real constant v. As mentioned before, spatially flat space is assumed, although 
there exist curved models in which the dynamics is also solvable by means of elementary functions [13, [38| . We define 
a scale factor x, so that it is equal to unity at time to (say, today), and accordingly, the scaling law for the (dust) 
matter density 

a(t) = a x(t), p(t) = p m0 x(ty 3 . (5) 

The main equation is now 



x'{ty A SllGpmO 



x(t) 2 3 

Next, we introduce the density parameters fl 



x{t)- A +upi n0 x{t)-\ (6) 



pm0 = ^G Qm0 ' u P™o= H fa<i, A = 3ff 2 O A , (7) 
which requires a change of the time variable 

x(t) = y(H Q t) = y(s). (8) 
so that the equation is simplified to a well known form 

4^ = "a + n m ,oy(s)- 3 + n q y(s)- G , (9) 
y { s y 

with the condition binding the present values of densities 

i = q a + f2 TO)0 + n g . (10) 

Note, that in our notation these O's are all constants. Finally, introducing a new dependent variable u 



y(s) = tfl(sj (11) 

the right-hand side becomes a polynomial 

iu'(s) 2 = n A u{ s ) 2 + n m , u{s) + a q . (12) 

which can immediately be solved by the use of trigonometric functions only as 

, -2Cl m0 V^A + (fiA + - 4«a^ 9 ) cosh[3v^A(a ± s)] + (Qa - n 2 n + 4n A Q q ) sinh[3\/ni(a ± s)] ^ n ^ 
«(-) = > ^ •■ , (13) 

where a is the integration constant. Obviously such a general solution is not very clear, the more so because a is in 
general a complex number. That is why we simplify the main equation a bit more with another set of substitutions 
(provided that 0) 

u(s) = w(3^\Cl A \s) = w(t), f2q = g|0 A |, f2 ro o = A*|^a|, (14) 
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and obtain 



w'(t) = sign(A)it;(r)^ + (j,w(t) + q, 



(15) 



or, if A = 0, 



w'(t) = tt mfi w(T) + 1 - Sl m , , 



(16) 



upon using (|10[) . Depending on the sign of A, we have three main groups of solutions, which we proceed to describe 
in more detail. 

III. A > 

The curves described by equation (| 1 5|) in the plane (w,w') are either two hyperbolae or straight lines, as shown in 
Figure |T]). If A = fi 2 — 4q is negative, we have the hyperbolae situated above and below eo = — (jl/2 - they are really 
the same solution, only with the sign of time interchanged. The main equation always admits such a symmetry, but 
as we want a universe which is expanding at this moment, we choose the solution with w'(0) > 0. It is given by 



w{t) 



-fjb + (2 + fi) cosh(r) + 2 v /l + g + /xsinh(r) 



(17) 



and shown in figure (|2j). 




FIG. 1: The possible evolutional paths of the system when A > 0, depending on the sign of A — ~ 



intersected by the w = 1 line are admissible in a given set of the parameters, eo 



fj,' — iq. Only the paths 
At/2, ei, 2 = (-A»=F VA)/2. 



When A > 0, we have the other two hyperbolae, but as Q m .o and A are both positive, eo is negative, and so only 
the right branch can intersect with the w = 1 line. If e 2 = (/U + \/A)/2 > 0, a bounce without singularity appears. 
The formula is the same as in (|17|) . and is shown in figure ([3]). 
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FIG. 2: A > and A < solution. 



w 




FIG. 3: A > and A > solution. 



This is in fact the generic setup, as A is usually assumed to be positive, and the p 2 contribution negative. This 
means that q is negative, making A positive. For this reason we consider bounce the generic feature here. 

If A = 0, eo becomes a double root of the polynomial in the right-hand side of equation (fT5)) and accordingly a 
possible static solution. However, it is negative in our physical setup, so the only possibility are the straight-line 
solutions. They are different in that they asymptotically evolve from the static solution in the infinite past, as can be 
seen in figure Formula (fT7|) still holds. 
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FIG. 4: A > and A = solution. 



IV. A < 



The situation is much simpler now, as shown in figure ([5]) — the only possible curve is a circle when A = /i 2 — Aq 
is positive. The condition (fTO]) means that the circle must intersect the w — 1 line, unless A — 0, which corresponds 
to a stable, static solution. If e\ = (—fi — VA)/2 is positive, the evolution is a singularity free oscillation. Obviously 
only one formula is needed to describe this case 



1 



w{t) = - —fi + (2 + fi) cos(t) + 2yJ-l - q - Jt sin(r) 



(18) 



and a typical example is shown in figure ([5]). 



V. A = 



Finally, the simplest case with zero cosmological constant, where going back to equation (p~6|) . This time the solution 
is even simpler - polynomial in time 

9 

u(s) = 1 + 3s + -fl m . s 2 , (19) 

where, as before, we determined the integration constants by assuming u(0) — 1 and u'(0) > 0. This is also a bounce, 
with a minimum size of 

U min = 1 - 7^—, (20) 

which is only positive if f2 mi o > 1, or equivalently fl q < 0, corresponding to a negative p 2 term in the Friedmann 
equation. Qualitatively the behaviour here is the same as that in figure ([3]), only the expansion is not exponentially 
fast. 



VI. CONCLUSIONS 



We have studied the FRW models with an additional (1 + z) 6 term in the Friedmann equation. Astronomical 
measurements based on the relation H{z) and general dynamics of such models ignore the number and the particular 
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FIG. 6: A < and A > solution. 



type of contributions - as long as they scale the same way. Among these, we described: brane theory, non-Riemannian 
modifications and loop quantum cosmology. Their combined effect, as observed through the Hubble's relation, could 
in theory be null, due to the indefiniteness of sign. Which is not to say that they would be indistinguishable (or null) 
when other observations come into account. 

Assuming the non-zero term a~ 6 , it is still possible to find exact solutions, and classify all possible evolutional 
paths. It was shown that exact solutions can be expressed in terms of elementary functions and depending on the 
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sign of A we have three main families (with the exact formulas provided in the respective sections). 

For A = (section V) we obtain a simple algebraic solution which is a bounce, and for a negative p 2 contribution 
the minimal value of the scale factor is positive. Qualitatively it is the same as the one presented in Figure [3] 

When A < (section IV), we have a periodic behaviour (Figure O of the scale factor, but as above, the minima 
are not necessarily positive. 

Finally, taking A > (section III), we obtain a universe which expand exponentially, but which could have followed 
three different paths in the past. It could have bounced (Figure [3]), it could have been expanding all the time from 
minus infinity (Figure ^ or asymptotically from the static solution (Figure 3]). Of course, the last two scenarios, 
would have to be considered as singular and expanding from the a = singularity. 

The generic feature seems to be the replacement of the initial singularity by a bounce. However, the effects might 
influence also the late evolution if one considers a phantom fluid, and are a possibility of avoiding the so called Big 
Rip singularity (39j . 

As mentioned before, similar exact solutions of the Friedmann equation were found for brane scenarios by Setare 
[UlUI]. Such solutions are among the simplest possible, where the left hand side of the equation is a polynomial of 
degree 2, or can be reduced to such a polynomial. Usually this can be done by either introducing a new "conformal" 
time, or choosing a new scale factor which some power of a. For polynomials of greater degrees such reduction is, in 
general, impossible and elliptic functions have to be used ■ 
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